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I. Introduction 


This report explains how a small digital computer, the Time/Data 100, 
that is specialized in the analysis of time-series data, is used to analyze 
magnetic field data. 

First, the computer will be described with reference to its operation 
and some of the algorithms it uses. Second, the organization of the data 
is detailed. Third, the results of the analysis of the data are examined. 
And finally, the Time/Data 100's method la compared with the Blaclonan- 
TuV.ey method of the calculation of the power spectrum. 
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II. Introduction to the Time/Data 100 Computer and to its Algorithms. 

The Time/Data 100 is a small digital computer that can process certain 
wired-in algorithms. These algorithms can be implemented by using the push- 
button controls on the console. 

Input data can be analog or digital and put into one of two 1001 word 
memory sections, each word having 7 bits plus a sign bit. The output data 
can also be analog or digital and is taken from one of two 1001 word memory 
sections in which each word consists of 17 bits plus a sign bit. Nine track 
IBM-360 compatible tapes are used for digital input and output. The input 
tapes are organized into 1011 byte records where the first 10 bytes are 
'header' data. The outp-c tapes also have 1011 byte records also with 
the first 10 bytes as 'header' data. 

The contents of one or two of the memory sections can be displayed on 
the cathode ray tube. All of the 1001 points can be shown or, with the use 
of the EXPANDED button, any 127 points of the memory can be displayed. 

The value of any word can also be displayed on the binary lamps. For 
example, if the contents of the memory section in INPUT A is being displayed 
on the lower beam of the tube, then the WORD SELECT switch is set to DISPLAY 
LOWER and the appropriate word is set on the digit switch. The value of that 
word will then appear on the binary lamps. It can also be seen on the tube 
by pushing the MARKER button. This will brighten 127 points, the first point 
of which is the selected word. If the EXPANDED button is pushed, these 127 
points will appear on the tube. 

When any of the 18 bit output memories are displayed on the tube, the 
AMPLITUDE SCALING switch must be used because only 7 bits plus the sign bit 
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of any word can be displayed. When the switch is set, the 7 hits shown on 
the tube are the bit number selected pi.ua the next six higher valued bits. 

For example, if the switch is set on 3, then the bits displayed are 2 s through 
2 9 plus the sign bit. 

There is also a plotter that will duplicate on paper any one of the 
memories that is displayed on the tube. Again, it can only display 7 bits 
plus the sign, so the AMPLITUDE SCALING switch must be used to select the 
dynamic range to be plotted. 

Connected to this processor is the system controller. It controls the 
data flow and enables the operator to process data continuously. It also 
has the algorithm and memory select switches that the console has. It shows 
the ’header 1 data and has a switch for manual or automatic operation. When 
the computer is in manual, the operator must execute each step himself and 
also select the correct memory section that the algorithm needs. For auto- 
matic processing, the computer will automatically select the correct memory 
section and will input, process and output the data. The processing will 
continue automatically until the MASTER CLEAR button is pushed, the limits 
of the algorithm are met, or until the number of times the algorithm was 
processed Is equal to the number in the REPEAT digit switch. This switch 
can be set from 00 to 99, 00 meaning Indefinite operation. 

There Is also a SEARCH switch that allows the operator to find a record 
on the input tape that is greater than or equal to the record number set in 
the RECORD digit switches. 

The three algorithms we are using are the Direct Fourier Transform, the 
Auto-Spectral Density and the Ensemble Average. 
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The Time/Data 100 uses this form of the Direct Fourier Transform for a 
discrete time series and finite time interval. The formula is: 

+N 

5 ( K f 0 ) = HI j- (nto) cos m k k ( n t o) 

h= -f\| 

Q(Kf„)- -[ L K h to) sin att f 0 (n t o)l 

L h: -N J 

where: S(Kfo) and Q(Kfo) are the real and imaginary components of the D.F.T. 

respectively. 

O S « £ IOOO 

+ 1 = number o'f dci+a, points 
i nq rat e 

'So - f n / 1 goo ^recjuehcy *=■ Yx \* o 



The final form of the algorithm Is the Rapid Fourier Transform that 
uses various time and frequency folds to shorten the calculations. The 
algorithm takes one Input frame and outputs two frames of data containing 
the real and Imaginary components of the D.F.T. 

The Auto-Spectral Density algorithm uses the following formula: 


where 0 £ K £ 1000. 


PfKfo) - | S(Kfg ) | 8 + | Q(Kfo)| 8 
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The algorithm takes 7 bits plus the sign, specified by the AMPLITUDE 
SCALING switch, of each word of the real and Imaginary components of the 
D.F.T. It squares and adds the two and puts It Into one of the output memory 
sections. 

The Ensemble Averaging Algorithm then adds together each successive 
output frame and puts It Into an output memory section. 

There Is, however, one algorithm that uses all three. It Is the 
Averaged Auto-Spectrum algorithm. It takes the D.F.T. of an Input frame, 
scales the two output frames, takes the sum of the squares and adds each 
successive output frame to form the averaged auto- spectrum. 

This algorithm can be run in a single cycle or continuous mode. When 
done continuously, It can accumulate at least 512 records and will stop 
automatically If there is an Impending overflow In one of the records. 

To test the algorithm, I made up a tape containing records of this 
function: 

Kt) r s \ l0sin ) 1-5 sin ( ^ + + 

Assuming the sampling rate is 1 second, then the corresponding Nyquist 
frequency is f N * -^-Hz, At ■ 1 sec. 

If the algorithm does work, then I should find three spikes at 0.083 Hz, 

0.166 Hz, and 0.333 Hz corresponding to the three sine functions. 

First, I used the D.F.T. algorithm to determine the correct amplitude 
scaling so that the power spectral density algorithm will use the correct 
values. An amplitude scaling of 7 was sufficient for the power spectral 
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density algorithm to take the highest value of the D.F.T. of the input funct- 
ion. Second, I ran the machine on automatic and accumulated 45 records. 
Third, I plotted the resulting averaged power spectrum and found out that 
there were three spikes at the expected frequencies. 
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III. The data and Its preparation for analysis 

The data we are using is from the OGO-E rubidium magnetometer experiment 
that measures the magnitude of the intensity of the magnetic field. The 
sampling rate is one point per 144 milliseconds. 

In order to do the power spectral analysis we are assuming that the 
data satisfies the ERGODIC HYPOTHESIS. First, it is a stationary, random 
process or the ensemble averages of the data don't vary with time. Second, 
each record is statistically equivalent to every other record. And, third, 
the ensemble averages of many records can be replaced by one typical time 
average. 

The data was sampled at three different rates: 

1. Every point or one point per 144 milliseconds. The corresponding 
Nyqulst frequency is 3.47 Hz. 

2. Every seventh point or one point per 1.008 seconds. The Nyqulst 
frequency is 0.496 Hz. 

3. Every fourteenth point or cne point per 2.016 seconds. The Nyqulst 
frequency is 0.248 Hz. 

I wrote a program to read the data tapes, sample them at these Intervals 
and build 1001 point records for input to the Time/Data 100. Each record is 
1011 bytes long where the first 10 bytes of header information contain the 
time of the first data point in each record, the record number, the ident- 
ification number and the scale factor. The average of each record vas sub- 
tracted to get rid of the DC component. 

The output tape containing all of the data sampled at the three rates 
was then used as input to the Time /Data 100. 
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To select the correct amplitude scaling for the power spectrum, the 
D.F.T. of each record had to be checked so that the most significant bits 
were used. Without doing this first, some significant bits, especially in 
the lower frequencies, would be lost and erroneous results would be obtained. 
After the amplitude scaling setting and been determined, the spectral 
averaging was done. 
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IV. Time/Data 100 Results 

Our examples are taken from October 26, 1968 from 6 H 0 M 15 S to 14^8^0^. 
The first graph Is a sample of the first record of the data and the second 
graph is its corresponding power spectral density. The curve is very rough 
but it can be smoothed by accumulating many record! The third graph is 
the accumulated power spectrum of the above time interval sampled at 144 
milliseconds. The resulting curve is much smoother than the power spectrum 
of one record. The fourth and fifth graphs show the time interval sampled 
at 1.008 seconds and 2.016 seconds respectively. These are not as smooth as 
the third graph, because of fewer accumulated records, but they still show 
the same general shape of the spectrum or the third graph. 
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V. Blackman-Tukey Method 

We are also comparing the Time/Data method of computing the averaged 
power spectrum with the Blackman-Tukey mecfcod. They express the power 
spectral density in terms of the autocovariance of the time-series. 

P(f)-- r“cCr) e - LiI ' fT ' d -, 

- *4 

w^e y 6 ; r 

c(t) =■ -f; m „ f J. it 

2 . 

We are using the following algorithm for our data: Given one record 

of the time-series f^, t = 0,1, ... ,1000; the total number of observations 
are 1001. 

1. Prewhiten the data. 

r St “ 0.(c ft-1 l * t i (OOO 

2. Calculate the autocovariance and subtract the mean. 



where r ■ 0,1,..., 100 ere the numbers of lege. 

3. Calculate the finite cosine series transform. 

V* = [Co + i E. + ^ io ° t06 
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4. Use of the ’‘Hanning" lag window is equivalent to the convolving of 
the computed power according to: 

u 0 - \ (V 0 + V t ) 

U r - i V r - 1 + j V r + J- Vr + 1 1 s r fi 99 

«ioe ' 2 V »» + 2 V ioo 

5. Correct for prewhitening and correction for the mean 


3b - I .I 0 COS t /30O 


a 



_J r \\Xr 

1.20 C04 fx/lOO J 

j 

- 1.20 cosE ( I " '/ bo®) ^ 



U-loo 


- 12 - 


VI. Conclusion 

We are presently debugging the program that prepares the data for the 
Blackman-Tukey method of power spectral analysis and will eventually be 
using both methods to analyze the data. 

Since the Auto-Spectral Averaging algorithm of the Time/Data 100 can 
process one record in one second, we can accumulate many records in a very 
short time. We can also selectively accumulate records and can compare 
their results in just a few minutes. So, the speed of computation and the 
immediate return of results are important factors in using the Time/Data 100 
for power spectral analysis. Therefore, with these two factors in mind and 
from the results we have already obtained, we feel confident that the 


Time/Data 100 computer can be successfully used to analyze our data. 
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Figure 3. Interface controller front panel 
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